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In this paper we consider a class of higher dimensional theories consisting of D-dimensional 
gravity coupled to a scalar dilaton and a form field propagating over a warped higher dimensional 
spacetime. Although in such theories it is possible to find exact cosmological de Sitter solutions, 
these suffer from instabilities that have been related to the runaway behavior of the moduli potential 
along the warped direction. In the simplest set-up, the models are characterized by two moduli: 
one related to the volume of the internal space, the other to the modulus of the warp factor. 
While the volume-modulus can be fixed by appropriately tuning the gauge field strength, curvature 
of the internal space, and cosmological constant, the same mechanism cannot work for the warp 
modulus. Here, we will present a stabilizing mechanism for the warp modulus and its mass in terms 
of quantum fluctuations from both moduli. We will show that, while quantum effects from the 
modulus associated to the warp modulus can only provide a stabilization mechanism of the mass 
scale in a restricted region of the parameter space, quantum effects from the volume modulus offer 
an efficient mechanism of stabilization. 

PACS numbers: 04.50.-h, 11.25.Mj, 98.80.Cq 



I. INTRODUCTION 



The recent discovery of dark energy demands a mecha- 
nism setting the cosmological constant to a value that is 
nonzero but hierarchically small compared to the Planck 
scale. At the same time, a lot of recent observational 
data, in particular those of the cosmic microwave back- 
ground, support the basic predictions of inflationary sce- 
narios. A theoretical framework that may be able to pro- 
vide a consistent description of the universe, undergoing 
inflation at early times and dominated by dark energy at 
the present day, is offered by string theory. In general, 
string theory requires the presence of extra dimensions 
that have to be stabilized at some appropriate scale to 
obtain a viable cosmological model. The lack of such a 
mechanism is often called the moduli stabilization prob- 
lem. The stabilization of moduli is deeply connected with 
the realization of above accelerating phases in the cosmic 
history, and plays an important role in the construction 
of higher-dimensional cosmological models. 

While the problem of moduli stabilization has been 
discussed at length for the case of Kaluza-Klein com- 
pactifications, its analysis for warped compactifications 
remains much less extensive I n fact, the case of 

warped compactifications is an interesting set-up to con- 
sider since, due to the warping and external directions, 
multiple scalar moduli may appear and mix in a non- 
trivial way, affecting, in principle, each other's dynam- 
ics. The present paper aims at discussing a rather sim- 
ple example of this sort. In particular, in our model, 



the corresponding exact solution of the warped de Sitter 
compactification in higher-dimensional gravity is known 
and the runaway behavior of the modulus implies the 
instability of the higher-dimensional solution. 

The class of models central to the present discussion 
consists of a warped higher-dimensional theories. Specif- 
ically, we will consider gravity coupled to a scalar dila- 
ton and a form field propagating on the background of a 
higher dimensional warped geometry of the form 

ds 2 = c 2A{ - v) [ds 2 (X) + ds 2 (Z) + dy 2 ] , (1) 

where X and Z are maximally symmetric manifolds and 
y is the direction of warping. Details will be given later, 
here we simply mention that an effective theory can be 
directly derived by compactifying the Z space, leaving, 
in the simplest construction, two unstabilized moduli: 
-0, related to the volume of the internal space Z, and 
A, related to the warp factor. As we will see, simulta- 
neous stabilization of both moduli may not be achieved 
by means of the same mechanism. For instance, appro- 
priately tuning the gauge field strength, curvature of the 
spherical internal space, and cosmological constant may 
help to achieve stabilization of the volume modulus ip 
but not of the warp factor 0, H[ ■ When adopting a sta- 
bilization mechanism based on quantum effects, even in 
the case in which vj is stabilized due to the presence of 
a gauge field strength, its quantum fluctuation as well as 
those of the modulus associated to the warp factor may 
both contribute to stabilize or destabilize the background 
geometry. 

At this point there are a few remarks to make. First of 
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all, after the stabilization of the volume modulus ip, quan- 
tum fluctuation of the volume modulus occurs around 
the minima that the effective potential attains in the ■in- 
direction. On the other hand, at tree level, the dynamics 
of A is controlled by a runaway type of potential, there- 
fore its quantum fluctuations should be analyzed with 
care. In the following we will use the background field 
method and tacitly assume that the system evolves adi- 
abatically. For this reason it is essential that the mod- 
ulus A evolves in a region where the potential allows it 
to 'slowly roll'. A second important point to remark is 
related to the value of the scalar potential after stabiliza- 
tion. In principle, once quantum corrections are included, 
the scalar potential of the system may occur at a posi- 
tive, vanishing or negative value, resulting in a de Sitter, 
Minkowski or anti dc Sitter geometry. In this case, we 
may expect that additional corrections to the potential, 
for example due to finite temperature effects, may pro- 
duce a further shift up-lifting its minima from anti-de 
Sitter to Minkowski or de Sitter, or, at very high tem- 
perature, pushing the system into a destabilized phase. 
The purpose of this paper is to analyze the possibility 
to achieve moduli and mass stabilization and explicit de 
Sitter solutions in the four-dimensional effective theories 
in the context of warped compactifications. 

The paper is organized as follows. In Sec. [Ill we will 
present the model in detail and construct the effective 
theory tuning the field strength to achieve stabilization of 
the volume modulus at the classical level. The main part 
of the paper is devoted to discuss how quantum effects 
from both moduli contribute to the effective potential at 
one-loop. We will adopt the background field method 
and path integrals to perform the computation and use a 
zeta function scheme to regulate the divergences. Specifi- 
cally, Sec. lIIII dcals with the contribution from the modu- 
lus A and analyze whether its quantum fluctuations may 
provide any stabilization. In fact, due to the runaway 
behavior of the potential, as we have already mentioned, 
the minima (if any) generated by quantum fluctuations 
must be in a region where the potential is flat enough 
and the perturbation theory can be trusted. This self- 
consistency condition will possibly limit the range of va- 
lidity of such a stabilization mechanism. After presenting 
the machinery we will perform the computation using an 
approach based on contour integral techniques similar to 
that described in Refs. (Related work is that of 

Refs. [IM1). 

This method is valid over the whole pa- 
rameter space and serves as a general way to compute the 
one-loop effective potential. In a restricted range of the 
parameter space a slightly simplified approach based on 
the Schwinger-De Witt approximation can be adopted. 
This method uses directly the small-i heat-kernel asymp- 
totics and it applies only in a small region of the parame- 
ter space. Details of this second approach will be given in 
Appendix [A] where the validity of the Schwinger-De Witt 
approximation will also be discussed. Results using both 
method are consistent when applied to the same region of 
the parameter space. (In Appendix [B] wc will show how 



finite temperature corrections may produce transitions 
between different minima uplifting the vacuum. These 
effects are studied by means of the standard Matsubara 
formalism. Clearly, this is only appropriate if the evo- 
lution of the system proceeds adiabatically, assumption 
tacitly made.) We will show that quantum corrections 
from the modulus A can stabilize its mass scale and at 
the same time provide either a de Sitter, Minkowski or 
anti de Sitter minimum. Unfortunately, the region of the 
parameter space for the moduli-stabilization consistent 
with the semi-classical approximation is only marginal. 
In Sec. II VI we present the warp modulus stabilization by 
computing the contribution which comes from quantum 
fluctuation of the volume modulus again using an ap- 
proach based on contour integrals. Quantum corrections 
from the volume modulus may provide an efficient frame- 
work for stabilization. This seems rather natural, since 
after flux compactification the size of the internal space 
generically becomes of order of the Planck length. In 
this case, it is not possible to ignore quantum fluctua- 
tions of the volume modulus, even though the volume is 
already stabilized. These contributions to the one-loop 
effective potential may stabilize the warped direction and 
naturally realize a de Sitter, Minkowski or anti de Sitter 
minimum depending on the values of the parameters and 
of the renormalization scale. Our conclusions close the 
paper. 



II. COMPACTIFICATIONS WITH MATTER 
FIELDS 

A. Compactifications with field strength 

The class of models that will be considered in this work 
is described, in the Einstein frame, by the following ac- 
tion 



S 



1 

2^2 



{i? - 2c- Q * /(p - 1) A} *l-^d(/)A *d<j) 



2-p\ 



e °"i>F A *F 



(2) 



where k 2 is the D-dimensional gravitational constant, * 
is the Hodge operator in Z?-dimcnsions, is a scalar field, 
F is a p-form field strength, and A and a are constants. 

The D-dimensional background geometry is assumed 
to have the following warped structure 



ds 2 = c 2A(y) [q^{X)dx»dx v + dy 2 
+ lab (Z)dz a dz b ] , 



(3) 



where is the n-dimensional metric over the manifold 
X and it depends only on the n-dimensional coordinates 



parameterizing X, while 7 a (, is the p(= D 



1) 



dimensional metric over the manifold Z that depends 
only on the p-dimensional coordinates z a parametrizing 
Z. The p-form field strength is taken to be proportional 
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to the volume form of Z, 

n(Z) = y/jdz 1 A • • ■ A dz p . 

that is 

F = fQ(Z), 



(4) 



with / constant and 7 denoting the determinant of the 
metric ~j a b- The choice (HJ guarantees that both the 
Bianchi identities and the equation of motion for the 
gauge field are automatically satisfied. 

The equations of motion for the scalar field and the 
Einstein equations can be obtained directly from @. As- 
suming the scalar field 4> to be proportional to the warp 
factor A(y), 



<P=-(p-l)A(y), 



a 



(5) 



and using ([3]) and (01, the equation for 4> can be written 
as 



A" + (D - 2) (A') 



p-1 



0. 



(6) 



A f 2 

where A = — ^Aj- + i-. Using the metric ansatz ([3]), 
Einstein equation can be expressed as 



R^(X)-/3Aq^(X) = 0, 



(7a) 



A" 



2{p-lf 



(P-1) 



(D-l)a>M -2^A = 0, (7b) 



R ab (Z)- (/3A+^-) 7 a 6 (Z) = : 



(7c) 



where ' denotes the ordinary derivative with respect to 
the coordinate y. i? MJ/ (X) and i? a (,(Z) are the Ricci ten- 
sors of the metrics and 7 a fc, respectively, and the 
constant (3 is defined by 



P = 



a 2 2(p-l) 



(8) 



p-1 D-2 
Eqs.© and (|7b| can be simultaneously solved as 

^(y)=^(y-2/o) , (9) 

with yo constant and £ given by 



± a ' 



A 



(p-l)(U-2) 



(10) 



Notice that the above solution is only the compatible 
with the condition a 2 ^ 2(p — l) 2 /(D — 2). Choosing a 
such that (3 > ensures that both X and Z are positively 
curved, as it is clear from an inspection of Eq. (|7c|). This 
corresponds to taking 



a > 



D-2 



a < 



L> - 2 



(p-1). (11) 



In this case, the field equations lead to the following so- 
lution for the D-dimensional metric 

ds 2 = c 2 ^y-y») [-dt 2 + <? m 8 pa dx p dx° + dy 2 

+ lab {Z)dz a dz b ] , (12) 

where the Hubble parameter H is given by 



H 2 



/?A 



(13) 



In the limit of H — > 0, we see that the Ricci tensor of 
the Z space leads to i? a b(Z) — > ^/ 2 7 a b(Z) . As the inter- 
nal space is essentially supported by the field strength, 
its properties remain unchanged in the limit H — > 0. 
From Eq. ([7a)l . so long as (3 > 0, a positive cosmolog- 
ical constant is allowed for A < ^j^-f 2 , which clearly 
means that the upper bound on the cosmological con- 
stant is given by gauge flux. As the final remark, for 
(D, p, a) = (10, 3, -1) with A = 0, the action © is 
that of ten-dimensional type IIB supcrgravity. For this 
choice of parameters, j3 = and the only possible solu- 
tion is (n + l)-dimensional Minkowski spacetime. 



B. The effective theory with field strengths 

In order to construct explicitly the effective theory, we 
can express the Z?-dimensional metric as 

ds 2 =e 2A(v) [ WpQ (M)dv P dv Q 

+e 2 ^ 1 ab(Z)dz a dz b \ , (14) 

where ip{v) represents the modulus degree of freedom 
of the Z space. wpq(M) and 7 a fc(Z) are, respectively, 
the (n + l)-dimcnsional metric and the metric of the 
(D — n — l)-dimcnsional Einstein space Z. In contrast 
to the previous subsection, the y direction is regarded 
as an external direction in the present effective theory 
approach. 

In the following, the (n + l)-dimensional effective the- 
ory is obtained after compactifying the Z space. In 
the matter sector, we will fix the gauge field strength 
F as in Eq. ([4| and require the scalar fields to obey 
Eq. ©. Hence, the modulus associated with the metric, 
namely the internal modulus ip and the warp factor A, 
are only the dynamical degrees of freedom in the lower- 
dimensional effective theory. 

Using ([5]) and HH), in the (n + l)-dimensional 
action ([2]), we get 

S = ^2 f [{R(M)-V(A, V5)}*m1m 



1 - - 1 C 2 - 

-dA A *mgL4 - - dA A *-M.dij) 

2 2 Jcic z 



1 



-dip A 



(15) 
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where i? a f,(Z) = A-f a b(Z) with A constant and the met- 
ric wpq(M) has been conformally transformed according 
to w PQ {M) =_e- 2 ^ D '^ A +P^/^-^wp Q {M). The Hodge 
operator on M space is defined as and k is given by 
K = V~ x I 2 k, with the volume of the internal space Z given 
by 



V= / * z l 



(16) 



The Ricci scalar computed with respect to the met- 
ric wpqiM) is R(M). In obtaining Eq. (fl~5|). we have 
dropped the surface terms coming from A M A A M 0, 
where A M is the Laplace operator constructed from the 
metric u>pq(M). The potential V (A, ip) is given by 

V (A,ip) = U(A)W(i>) , (17) 

where 



U(A) = exp 
W§) =2Aexp 



2(D - 2) A 



(n - l)VcT 
2pip 



(n - l)Vc3 



exp 



2npip 



(n - 1) V / C3 



pX exp 



2(D - 2)V> 
(n - ^v 7 ^" 



(18a) 
(18b) 



The fields ^4, ^ have been rescaled according to 



III. STABILIZATION OF THE MODULUS MASS 



A = JcyA , 1p = yjcii) 



(19) 



with the constants Cj (i = 1, 2, 3) defined by 



ci = 2 



C2 



Ji- 1 

+2 

4(D-2)p 
n- 1 



(D-2)- 2(D - 1 



n-l + -j(p-l) 



, n - 1 
c 3 = 2p [ + 1 



(D-2) 
(p - 1) + 2p(D - 1) , (20a) 
(20b) 
(20c) 



The absence of a stabilization mechanism for the modulus 
associated to the warp factor is clear from the form of the 
potential in Eq. (fT7|) . In the A-direction the potential de- 
cays exponentially causing the modulus A to suffer from 
a runaway behavior and the warped direction to expand 
forever. This runaway behavior is consistent with the ap- 
pearance of an unstable scalar mode found in Ref. 0, Q ■ 
In the present set-up, classically, the warped direction 
cannot be stabilized. On the other hand the vacuum ex- 
pectation value of ip can be fixed by appropriately tuning 
the gauge field (see Fig. QJ. There is certainly a local 
minimum in the ip directions. The potential energy at 
the minimum is equivalent to the (n + l)-dimensional 
cosmological constant. Since the moduli potential en- 
ergy eventually turns out to be positive or negative, the 
(n+l)-dimcnsional background geometry becomes dS rl +i 
or AdS ra +i spacetime. In the following, we will assume 
that the modulus "0 is fixed at ij) = t])o by tuning the 
gauge field flux. This does not affect the dynamics of the 
other modulus A, whose stabilization will be considered 
in the next section. 



In this section, we will discuss the possibility of stabi- 
lizing the mass scale of the modulus degree of freedom 
associated with the warp factor in the lower-dimensional 
effective theory incorporating the solution described in 
Sec. UTAl 

The goal is to determine whether quantum fluctua- 
tions of the modulus A can provide a reasonable stabi- 
lization mechanism of its mass scale. Here, we will adopt 
the background field method and the path integral ap- 
proach to compute the effective potential at one-loop for 
the moduli-field and deal with the divergences using zeta- 
function regularization. 

Let us consider the (n+ l)-dimcnsional scalar sector of 
the action (fl~5l) 



1 



M 



1 



: dA A *m<L4 



+V (A, ij ) * e 1m] , (21) 



and expand the field A around its classical vacuum ex- 
pectation value, Aq, 



A(v) = A + a , 



(22) 



with a representing the quantum fluctuation. Expanding 
the action up to second order 



S = ^2 I [{R(M)-V{a, A , 0o)} 



'M 

— —da A * M c?a 



(23) 



where linear terms in a(v) have disappeared owing 
to the classical equations of motion and the potential 
V (a, A 0l ipo) is given by 



V (a, Ao, fa) = U(A )W(fa) 
+0(a 3 ). 



2(D-2) 2 2 
[n — 1) ci 



(24) 
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FIG. 1. The figure illustrates the behavior of the potential W(ip). The left panel displays the potential for vanishing gauge 
field, / = 0, and for several values of the constant A (with A normalized to unity). The right panel shows the potential for 
several choices of / (with the values set as indicated in the figure). Dimensionality parameters are chosen as follows: D — 10, 
p — 6 and n = 3. 



After varying the action with respect to a, we obtain the 
field equations for the fluctuation 



(Am - Ml) <* = 
where Ml is given by 

r2 _2(D-2f 



(25) 



Mi 



U(A )W(^o) 



a (n - l) 2 Cl 
Using path integrals we can express the amplitude as 



(26) 



Z = J V[a}exp(iI A [ip , A]), 



(27) 



where T>[a] is a measure on the functional space of scalar 
fields a(v), and IaI^o, A] is given by ([21]) . At one- loop, it 
is sufficient to compute the above path integral with the 
action expanded up to second order around its classical 
background value, 

I A [4> , A(v)] = I c [fa, Aq] + I q [fa, A , a] + 

+0 (a 3 ) , (28) 

where A(v) is given by (f2"2"j) and linear terms in a have 
disappeared due to the classical equations of motion. Us- 
ing the above expression, the path integral (|27j) becomes 



lnZ = U c [ipo, Aq] + 

+ ln ( f V [a] exp (il q [^ , A , a] H . (29) 



The above integral is ill-defined because the operators 
in Eq. (|29p arc unbounded from below in the dS n _|_i 
spacetime with Lorentz signature. In order to correct 
this pathology, we proceed in the usual way and by per- 
forming a Wick rotation re-express ([2l?|) in the Euclidean 
form, 

\nZ= -I cE [tpo, A Q ] + 

+ ln( f 2?[a]exp(-/ qE [V;o, A), a})}, (30) 



where J q E is the Euclidean action expressed by 

IqE^o, Ao, a] = — g / a [-d *m d + Ml *m 1m] ^S 1 ) 
4K Jm 

Here we have integrated by parts over the kinetic term. 

The one-loop quantum effective potential V q is defined 
according to the relation 



exp 



- / Vq(Vo, A ) 



■M M 



M 



= / V [a] exp (-IqE^o, A , a}) 



[det yT 2 ( A 



M 



Ml)] 



(32) 



where A M denotes the Laplace operator on (n + 1)- 
dimensional dc Sitter spacetime, and fi is a normalization 
constant with dimension of mass. Defining 



M 



exp - / V q (ip , A ) * m 1 m = exp [-fl V oiV q (ipo, A )] , 

(33) 

with fl vo \ being the volume of (n + l)-dimcnsional de 
Sitter spacetime, we obtain the following expression 

Vcft^o, A ) = V (tp , A ) 

lndet [ M - 2 (A (M) - Ml)] , (34) 

Zi £ V ol 

where the above functional determinant has to be evalu- 
ated on dS ra +i. 

A natural way to proceed is to use zeta regularization 
techniques. Defining the following generalized zeta func- 
tion 



(35) 



where A are the eigenvalues of the Laplacian on dS„+i, 
the effective potential ([34]) can be expressed as 



VeffW'O, A ) = U (^o, Aq) 

~^[Ca'(0) + 2C a (0)ln( M 6)], (36) 
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where b is the radius of a (n + l)-dimensional sphere 
S n+1 . The task is then to find the analytically continued 
values of the zeta function and its derivative, Co(0) and 
C a '(0). The one-loop effective potential can be computed 
in a variety of ways. The most advantageous one is to 
use contour integral techniques, which will be done in 
the reminder of this section. However, to see the overall 
feature of the effective potential, the simplest way would 
be the 'Schwingcr-De Witt' approximation, which will be 
performed in Appendix ([A"]). 

The (n + l)-dimensional de Sitter geometry, dS n+ i, is 
a (n + l)-dimensional manifold with constant curvature 
and has a unique Euclidean section S™ +1 with a radius b. 
We call the eigenvalues of the Laplacian on this spacetime 
A (I) and their degeneracy d{£). These are explicitly given 
by El 



where we have defined 



Dn ( L - - 



2L-1 
(2N - 1) 



1 

L -2 



L 



L -2 ] 



N -2 



(41a) 
(41b) 



Using the residue theorem, we can replace the infinite 
mode sum over L by complex integration, obtaining 



d{£) 



{21 + n)(£ + n- 1)! 



\(£) = £(£ + n). (37) 



Using the generalized zeta function Eq. (|35p which can 
be explicitly re-expressed as 



b 2 



Mi 



(38) 



Z±{B) = --( — 



2 s 



B N x 



x / dz ta,n(B N irz)D n (B N z) (z 2 =p l) 



,(42) 



where the contour Ci in the complex plane is showed in 
Fig.fJJ and B^ is defined by 



the effective potential is given by Eq. (|36|) . We will eval- 
uate the analytically continued values of the zeta func- 
tion (|38[) at s = referring to the method employed in 
Refs. (ail. 

We perform the analytic continuation of the general- 
ized zeta function to s = in the case of n being an odd 
positive integer, since the value we are interested in is 
?i = 3. Then, 



~) -(bM a ) 2 



(43) 



(For a positive B^, -(N - \) < bM a < N - |). For 
clarity, we will consider the two cases separately, 



Ca(*)=E 



e=o 



(2£ + n)(£ + n-l)\ 
ri£\ 



• + n) 



IP 



Mi 



(39) 



Defining N = (n + l)/2 and using it as running variable 
L = £ + N, we rewrite the above expression as 



Us) = J2 D « 



L=N 



1 



A„ L 



b 2 



Mi 



(40) 



, , Z + (s) if B* >0, 
- i zS{s) if Bl < 0. 



(44) 



Let us consider Z+(s) first. In order to avoid the branch 
points z — ±1, we may proceed by deforming the contour 
Ci into C2 as indicated in Fig. [2] (left panel), and express 



as 



Z+(s) = 



b_ 
B N 

i 

2 C 



2 s 



1 



dx D n (iB^x) (x 2 + l) tanh(i?N7i"^) 



dx tax\{B^Tr{x — ie)} D n (B^x) (l — x 2 ) 

r 1 

-e _I7rs / dx tan{B N n(x + ie)}D n (B N x)(l-x 2 ) 
2 Jo 



(45) 
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where D n (iB^x) defines the following polynomial with coefficients r^k 



D n (iB^x)=i{-lf^ J ^^_ 



x) 2 +[N-^ 



N-l 



k=0 



K-lf- 1 £ r Nk (B N x) 2k+1 . 



(46) 



The first term in Eq. ([45)) comes from the integral along the imaginary axis. The second and third terms in Eq. ([45]) 
arc the contributions from the contour along the cut on the real axis. Using in the first term of (|45p the following 
relation 



tanh(i?Ni'2;) = 1 



(47) 



we arrive at 



Z+{s) 



b 



N-l 



DsM-Lt-i)^- 1 £ r Nk (B N ) 2fc+2 r(fc + i)r( s - *- 1) 

f°° - s 2 

i -Bn cos(7rs) / dx D n {iB^x) (x 2 + \) 
Jo 



g2_B N 7T.X' _|_ ^ 



i?N sin(7rs) / dxD n {B^x) (l — x 2 ) s ta,n(B^nx 



(48) 



Next, we consider the function Z-(s). This time, the branch points in the integrand are on the imaginary axis at 
z = ±i. Therefore we deform the contour as indicated in Fig. [5] (right panel) and obtain 



Z_(s) = i 



\B N \ / dx D n (i\B N \x)tanh(\B N \nx) (l - x 2 )~ 




N-l 



{-irr{-s + l)Y,r Np {\BK\? P+2inS P 



p=0 



C0S(7TS) + 



r(p + i) 



2i|Bn| < cos(7rs) / dxD n (i\B^\x) 



x> - 1 



g2|B N |ira; _|_ ^ 



r(-p) 



3 ^^(^Nk)- e2|BN| ^ + 1 




(49) 



The above expressions, (|48|) and (|49|) . can be easily ex- 
panded to get the analytically continued values Ca(0) and 

C(o). 

In Figs. [3]|3l for n = 3, the behavior of V e g is numeri- 
cally illustrated as a function of M%, with three param- 
eters /z (or dimensionless ^6), a and b. In the left panel 
of Fig. [3j the effective potential is shown for various fib 
while fixing b = 0.4 and a = 1.0, and in the right panel it 
is shown for various a while fixing b = 1.0 and fib = 1.0. 
On the other hand, in Fig. 21 it is shown for various b 
while fixing a = 1.0, fib = 10. For a decreasing fib with 
fixed other parameters an AdS vacuum is lifted to dc Sit- 
ter or Minkowski one. If a is below a critical value for 
a given set of other parameters, it is not possible to find 
a vacuum. Finally, for an increasing b with fixed other 
parameters, the energy density of the de Sitter minimum 
decreases but the potential minimum eventually disap- 
pears before it becomes a Minkowski or AdS vacuum. 
Similarly, if there is an AdS vacuum, as increasing b it is 
lifted but eventually disappears 



I 

before it becomes a dc Sitter or Minkowski vacuum. 
The results of this subsection are confirmed by those ob- 
tained using the 'Schwingcr-De Witt' approximation as 
described in appendix |XJ 

In the next section, we will consider the contribution 
to the potential of the quantum fluctuations of the vol- 
ume modulus around the classical minimum ij) = ipQ de- 
termined by effects of the gauge flux and of the bulk 
cosmological constant. 



IV. QUANTUM CONTRIBUTION TO THE 
EFFECTIVE POTENTIAL FROM THE 
VOLUME-MODULUS 

In this section, we consider the case when the size of 
the internal space approaches the Planck length. In this 
case, quantum corrections can no longer be neglected. If 
the size of the internal space is larger than the Planck 
length, quantum effects can be analysed using the con- 
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FIG. 2. The left panel shows the deformation of the contour Ci used in (|42f> . replaced by the contour C2 that avoids the branch 
points at 2 = ±1. The right panel shows the deformation of the contour Ci used in (|49[) . replaced by the contour C3 running 
parallel to the imaginary axis. The points x^ are defined as: x^ — B^ 1 (N + k + 1/2) for (|45[> and x^ = B~^ 1 (N + k — 1) for 





0.5 Mi 



FIG. 4. The figure illustrates V c s for M| for changing 6 = 
0.8 — 1.2 from the top (red) and fixing a — 1.0 and /ib — 1.0. 



ventional loop expansion. In the opposite case, the loop 
expansion breaks down. Therefore, in the following, we 
assume that the radius of the extra dimensions is larger 
than the Planck length, which can provide a natural cut- 
off scale to the quantum field theory. 

Even if stabilized by flux, the volume modulus ip may 



still contribute to the dynamics of the moduli associated 
to the warp factor, A through the coupling of the quan- 
tum fluctuation of to A even if we ignore the quantum 
fluctuation of A itself. In this section, using the contour 
integral method, we will compute the quantum contribu- 
tion of the modulus at one-loop and discuss whether 
they can stabilize ip. As for the quantum corrections of 
A, we can expand 0(i>) around a neighborhood of the 
local minimum of the potential W{ip), 

0(0) = O + <p(v) , (50) 

where tJjq is fixed owing to the gauge flux (see Sec. II 
B). The (n + l)-dimcnsional action (p~5|) expanded up to 
quadratic order in ip(v) is 

/_ [{R(M)-V(A, O , ^)}*m1m 



5' 



2k 2 
1 



M 



-dtp A *ud(f 

where the potential V (A, 0o, ¥?) is given by 

V (A, 0o, <p) = U{A) [W (0 O ) + W 2 (0 O ) </? 2 
+t% 3 ) ■ 



(51) 



(52) 
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Note that linear terms disappear owing to the classical 
equation of motion and the second term explicitly de- 
notes the coupling of the quantum fluctuations of the 
volume modulus tp to the warp factor A. In the above 
expression, the functions Wi(^o) (* = 0, 2) are 



W (ip ) 
W 2 ($o) 



2Ae(p) 
2 



e(np) 



2A P 2 e( P )- 
-p\(D - 2) 2 e{D - 2)} 



p\e(D - 2) 
{npff 



(53a) 



e(np) 



where we have defined 



e{x) 



Varying the action with respect to (p gives 



A, 



o, 



with M 2 is expressed by the relation 

U(A) W 2 (4) • 



Ml 



(53b) 



(54) 



(55) 



(56) 



The calculation of the one-loop effective potential can 
be carried out using path integrals, and similar steps to 
those used in the previous section allow us to obtain 



V eS (A, Vo) 



U(A)W (^ ) + V q (A, Vo) 
U(A)W {i; ) + —J— In dot [p 



2fl 



vol 



2 ( A M 



M 2 V )} , (57) 



where Ajyj denotes the Laplace operator on (n + 1)- 
dimensional de Sitter spacetime, and fi is a normaliza- 
tion constant with dimension of mass. The (n + 1)- 
dimensional de Sitter geometry, dS rl+ i, is a (n + 1) 
dimensional manifold with constant curvature and has 
a unique Euclidean section S n+1 with a radius b. In the 
following we will evaluate the potential by analytically 
continuing the generalized zeta function 



e=o 



Mi 



(58) 



to s — > 0. The effective potential V c g(A, ipo) is then 
expressed as 

V cS (A, ij ) = U(A)W W ) 

-^-[C/(0)+2C v (0)ln( At 6)]. (59) 

Zi &vol 

We will refer to the method employed in Refs. [1, 
The contribution of the quantum correction played an 
important role to the effective potential. We find that 
the quantum effective potential has a terms proportional 
to M v . The procedure is the same as that employed in 
Sec. III. C, except for the replacement of M 2 — s- M^. 

As before, here we will focus on the case of n odd and 
integer. Using the residue theorem, and defining 



N~l) -(bM v ) 2 



(60) 



(B 2 ] is positive for -(iV - |) < bM v < N - |), we will 
consider the two cases separately, 



Z + (s) 
Z-(s) 



if B 2 > , 
if B^ < . 



(61) 



Following the same procedure as that in Sec. Ill C, we 
can finally reduce Z + (s) and Z_(s) to the same forms as 
Eqs. (|4"8"|) and (|4"9"|) . respectively, with the replacement of 
the definition of B 2 N as Eq. f|DJ) . 

The above expressions for the zeta functions can be 
directly used to obtain the one-loop effective potential. 
While explicit expressions can be obtained from for- 
mulae (|3H]) an d (SU), here we follow a more expedite ap- 
proach based on numerical approximation. Results are 
shown for the case of n odd that we set n = 3 and B 2 N 
positive. Figs. [5][7] illustrate the effect of the one-loop cor- 
rections from quantum fluctuations of the volume mod- 
ulus tp after flux stabilization. Depending on the value 
that the potential W(ip) attains at the minima, various 
possibilities can be realized. Fig. [5] shows a typical con- 
figuration that realizes a de Sitter minima, for positive 
W{^o) = 10 _1 . The right panel of Fig. [5] shows how the 
potential depends on the value of W(ifio) illustrating how, 
for increasing values of VF^o) the vacuum can be lifted 
from AdS to Minkowski or de Sitter. The quantum cor- 
rection basically lifts the potential up without changing 
the shape too much around the minimum, which allows 
to uplifts the AdS minimum and make it a metastable de 
Sitter ground state. Fig. [S] depicts the dependence of the 
potential on the parameter a, showing that a decrease 
in a tends to shift the minima towards larger values. 
Finally, for W{4>o) negative, again an AdS vacua is real- 
ized and increasing Wfyo) does not change the sign of 
the minima of the effective potential as long as Wfyo) 
remains negative (see Fig. [7]) and the potential tends to 
accumulate on the upper curve. For values of > W(?po) 
below a certain critical value, it is not possible to achieve 
any minima when quantum effects are included. 

The classical potential of A forces to decompactify the 
extra dimension while the combinations of matter and 
quantum correction produce a local minimum of the ef- 
fective potential. Hence, the scale of the internal space 
Z is stabilized by balancing the 1-loop correction, gauge 
field strength wrapped around the internal space and the 
curvature term of the internal space with the cosmolog- 
ical constant. If we can have a negative potential mini- 
mum for a choice of the parameters, a dS„+i spacetime 
evolves into a AdS„+i when the modulus settles down to 
the potential minimum. Associated with moduli stabi- 
lization, the change of spacetime geometry takes place. 



V. DISCUSSIONS 

In this paper, we have tackled the issue of the mod- 
uli stabilization in higher dimensional gravity. The class 
of models we have considered here presents two moduli. 
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FIG. 5. The figure in the left panel shows a typical configuration realizing a de Sitter minima after quantum stabilization. 
The small superposed figure represents the classical potential W(i[jo) for the volume modulus tp after flux stabilization. In the 
left panel we have set: A = 1.2, A = 0.547, / = 1.6 and a = 1. The right hand panel shows how the minima of the potential 
depends on the value of W(ipo)- The top green curve corresponds to W(ipo) = 1CP 1 realizing a de Sitter vacuum, while the 
bottom purple curve corresponds to W(ipo) = 10~ 2 realizing an anti de Sitter vacuum. The red dotted line corresponds to 
W(ipo) — 0.053 and realized a Minkowski vacuum. 




FIG. 6. The figure illustrates the dependence of the potential on the parameter a. In the left panel we have set the parameters 
A, A and / as in the previous figure in order to obtain, after flux stabilization, a positive minima, W(V>o) = 10~ . In the right 
panel we have reduced the flux to obtain W(ipo) = 10 -4 . In the first case (left panel), after quantum stabilization the minima 
is positive realizing a de sitter vacua, while in the right panel the minima is negative giving an anti de Sitter vacua. Decreasing 
the parameter a shifts the minima towards larger values, without changing the sign of the potential. 



One (-0) is related to the volume of the internal space, 
while the other (A) is related to the warped direction. 

These theories provide interesting cosmological toy- 
models owing to the fact that it is possible to realize 
explicit exact de Sitter solutions. But, unfortunately, 
it has turned out that these solutions suffer from an 
instability about the perturbations. In previous work 
(see Ref. [3, [||), the lower-dimensional effective theory 
has been derived, with the warped direction regarded 
as an external one and the warp factor as a modulus. 
This analysis suggested that the instability of the higher- 
dimensional solution was related to a runaway behavior 
of the warp factor in the lower-dimensional effective the- 
ory. Thus, the main goal of this paper was to discuss a 
consistent mechanism of stabilization for the warp factor 
and therefore a cure to the above instability. 

The example we have considered is simple enough, in 
the sense that only two moduli are included in the anal- 



ysis. While the volume modulus can be fixed by ap- 
propriately tuning the gauge flux, the same mechanism 
cannot work for the modulus associated to the warp fac- 
tor. Therefore, in the present paper, we have discussed 
whether quantum fluctuation from both moduli can lead 
to full stabilization. We have discussed this by using 
the background field method, path-integrals and zcta- 
function regularization, and showed that, while quantum 
effects from A can only provide a reasonable stabiliza- 
tion mechanism of modulus mass in a restricted region 
of the parameter space, quantum effects from ip offer an 
efficient solution to the moduli stabilization problem. In 
the presence of the 1-loop correction, the classical con- 
tributions of the curvature and flux terms compete with 
the quantum effects from the other ingredients; we have 
found the two to balance to yield a local minimum in 
this direction and showed that by tuning the choice of 
a and \i one can perturb the AdS vacua to produce dS 
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FIG. 7. In this plot we show how the minima of the potential in the A-direction shifts when the flux stabilizes the volume 
modulus to a negative value, generating an AdS vacua. The top curve refer to W(ipo) = — 1CP 4 while the bottom curve refers to 
W(tjjo) = —5 x 10 -3 . One may notice that for negative and decreasing values of W(i/jo) the minima of the potential accumulates 
around V m in ~ 0.0395 for the present choice of parameters. For values of W{tpo) below 5 x 10~ 3 the minima disappears. 



vacua with a tunable physical parameters. The vacua 
will clearly only be metastable, since all of the sources of 
energy we have introduced vanish or negative as A — > oo. 
In general, we have observed that depending on the val- 
ues of the parameters, either a de Sitter, Minkowski, or 
anti de Sitter vacua can be realized. 

If we replace the internal space Z =S D ~ n ~ 1 with 
QjS L/ , we can obtain another class of warped de Sitter 
compactifications. Since the solution has a similar form 
of the metric to [2lL l22j , more moduli which correspond 
to the scale of the S Ll appear after compactification. In 
these models, the moduli potential may contain a large 
number of metastable states, and the space of vacua is 
largely or completely connected, as in the cosmological 
landscape picture. Although a generic cosmological solu- 
tion in the landscape is complicated, a number of possible 
choices of vacua may be crucial for realizing the observed 
scale of the cosmological constant as proposed in [23| . 
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Appendix A: 'Schwinger-De Witt' approximation 

Here, we will provide a simpler way to compute 
the one-loop effective potential (|5rj|) directly using the 



Schwingcr-De Witt expansion for the heat-kernel. This 
approach is valid in the region of parameter space for 
which the value of M a is large enough. 

Using the Mellin transform, the zeta function can be 
expressed as 

Ca = ^ dtt^e-^+^Qit), (Al) 
r 0) Jo 

where H := b^ 1 is the Hubble scale of the de Sitter space 
and the function Q(t) is the heat- kernel defined as 

e(t) = ^ c - (A - ff2) *. (A2) 

A 

If the value of the mass M a is large enough, then the 
exponential in the integral above suppresses the contri- 
bution coming from the large-i part of the integration 
range, and a direct use of the small-i expansion is pos- 
sible. This procedure is analogous to the high tempera- 
ture expansion of the effective action. After rescaling the 
integral (|Alj) by t — > H~ 2 t, it is straightforward to real- 
ize that the exponential suppression becomes substantial 
when M 2 H~ 2 becomes large enough. Using ([26]) . it is 
straightforward to see that choosing a ~ 0(1) and tun- 
ing the gauge flux in such a way to obtain Wtyo) ~ O(l), 
a small hierarchy between the Hubble parameter H and 
the Planck mass (H ~ 10 Mpi) is sufficient to generate 
enough exponential suppression. 

In this region we may approximate the integrand in 
(|Aip by using the Schwingcr-Dc Witt expansion for 8(t) 

where the coefficients 9k are the heat-kernel coefficients 
[l5l [lr| . Explicit form for the coefficients can be found 
with little work and for the present case of de Sitter space 
with £,c — ja and n = 3, these are 
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6q — f^volj 

where f2 V oi is defined by 



h = -0 



vol 



f- 

V480 



= -n. 



/ 457 
V 40320 



if' 



vol 



(A4) 



(A5) 



A direct computation gives for the one-loop effective potential for D = 10 and n = 3 the following expression 



V eS (M a ) =V + V q (M a ) 
8a 2 + 25 



32a 2 



-Mi 



32tt 2 

, 2 



- {Mi + 2H J 



I'- 



ll 6 



71/2 + 2H 2 J 315 M 2 + 2ff 2 



(A6) 



where we have rescaled the various quantities according 
to 



V cff 
M a 



V eS k 4 

M a fi, 



H = Hk , 
p — fik . 



Eventual non-vanishing minima of the potential deter- 
mine the mass of the field a: 

= AV eS {a,p,H,M a ), (A7) 

0=-^AV en (a,fi,H,M a ), (A8) 

where we have normalized the potential according to 

AV e «(a,p,B,M a ) := V cS - 3H 2 . (A9) 

For a given set of (a, p, H), the solution for M a leads to 

V c « ~ 3k~ 4 H 2 . (A10) 

Hence, for H to be H < 1 and (M a /H) > 1, then 
if one keeps \x = O(l) fixed, the energy density at the 
minimum is much smaller than the Planck scale, which 
implies that the stabilization due to the quantum cor- 
rections is working consistently. In case of D = 10 and 
n — 3, the classical potential approaches a constant from 
above as a increases. For a tuned to be small but non- 
zero, the quantum correction no longer contribute to the 
effective potential. For modest values of a, we will find 
numerically that there is a solution M a ~ 12 for jl > 10. 

Approximate expressions for the minima of the po- 
tential can be found at leading order by expanding for 
M a 3> H . In this regime the minima is determined by 



Mi 



1 + ln 



7T i 



■8a 2 



(All) 



Assuming the renormalization scale to be of the same 
order as the mass, p, ~ M a , we find 



Mi 




a 2 [i 2 



(25 + 8a 2 ) 



(A12) 



Higher order corrections do not change the qualitative 
features of the above result. 

The value that the potential attains at the minima 
depends on the choice of the renormalization scale. Min- 
imizing Vmin as a function of /2 allows to find a Minkowski 
vacua (Vmin = 0) for 



Merit "|^( 25 + 8 " 2 ) = \ 



8tt 2 



(A13) 

where ^ m i n is the minimum value of the renormalization 
scale for which a minima with positive vacuum energy 
exists. An AdS minimum (Vrnin < 0) is found for values 
of p, in the range p > p cl -it, while a de Sitter minimum 
(Vmin > 0) is obtained for values of p lying in the range 
Mmin < p < pcrit an d the expansion rate is given by 



3H 2 



64tt 2 



1 



3tt 2 



2a 2 /i : 



■(25 + 8a 2 ) 



1 - 



a 2 /i 2 



(25 + 8a 2 



(A14) 



The above arguments, although apply in a specific region 
of the parameter space of the model (p and a) suggest 
that a can be stabilized by quantum effects. 

A more general computation of the one-loop effective 
potential valid in all regions of the parameter space was 
given in Sec. IIII1 which exhibits a behavior consistent 
with the results shown in this Appendix. The dependence 
of the potential on the energy scale suggests that the in- 
clusion of finite temperature effects may lift the minima 
of the potential. Of course, these effects are not directly 
related to the mechanism of stabilization discussed in this 
paper, and clearly a proper inclusion of thermodynamic 
effects requires care, particularly if time dependence is 
taken into account. However, in the approximation that 
the time evolution of the moduli fields is adiabatic, it 
is possible to give an estimate of these effects using the 
standard Matsubara formalism. The argument becomes 
simpler if the scale of the S™ is approximately constant, 
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i.e. if we assume the adiabatic expansion in the direction dius of the spatial section S™. The computation of the 
of S n after compactification of the (n + l)-dimcnsional potential at finite temperature carried out in Appendix 
theory over S 1 to S™ +1 =S 1 xS™ with H~ x being the ra- IBl gives for n = 3, 



elf 



25 + 8a 2 - 2 



1 



32a 2 

OO 

+ 4 E' 



32tt 2 

x(M 2 + ff 2 
ir 2 £ 2 



[Mi + H' 



3 1 ( Ml + H 2 

4 2 V 



'M 2 + iJ 2 



\ 



(A15) 



where \ = k 1 ^ 2 (2t:T) 2 and T is the temperature. Details 
along with high- and low-temperature approximation are 
obtained in Appendix [Bj Here, we show the typical be- 
havior of the potential in Fig.[8]where we have normalized 
its value by subtracting the vacuum energy contribution 
for M a = 0, which corresponds to the A — > oo limit. 



Appendix B: Finite temperature corrections 

In this Appendix wc present the computations of the 
finite temperature corrections to the effective potential. 
As mentioned in Appendix. [A] we assume that the time 

I 



r 



evolution of the modulus is adiabatic jjdlowing us to use 
the standard Matsubara formalism [201 . In this adia- 
batic regime, the scale of the S™, iT -1 , is assumed to 
be approximately constant. The same formalism of Ap- 
pendix [A] can be applied and the zeta function becomes 



(Bl) 



A. TV 



where A are the eigenvalues of the Laplace operator 
on the n-sphere, N is integer, and x is defined by 
X = (2irT) 2 . Mcllin-transforming the above expression 
and using the same heat-kernel scheme adopted in Ap- 
pendix \Xl it takes simple steps to arrive at 



as) 



(47r)2r( s )V x 



y/ X (Mj + H 2 ) 



[Ml 



K 



2ttA 



' M 2 + H 2 



X 



(B2) 



The important values £(0) and C'(0) can be computed in a straightforward manner from the above expression, leading, 
for n — 3. to 



eff 



25 + 8a 2 



32a 2 

OO 

+ 4 E' 



32tt 2 

X {Ml + H 2 
ir 2 l 2 



Ml + H 2 



K, 2ir£< 



3_i ln (Mj + m 

4 2 V P 



IM 2 + H 2 



X 



(B3) 



where we have used the definitions (|A7[) . and x = k 1 / 2 ^. The volume factor and the heat kernel coefficients arc given 
by 



2tt 
7* 



2fe-3 



(B4) 



Here is the volume of S 3 , and % is given by 70 — 1 and % — (k = 1, 2, 3) for an S 3 . 

Below we obtain the limiting behavior of JB3| assuming M a 3> H. At low temperature x — » 0, the modified Bessel 



functions decay exponentially as e 
the result of Appendix [A] 



■Jx , and the finite temperature corrections become small. Hence we can recover 
1 



25 + 8a' - a 
32a 2 a 



32tt 2 



4 2 



Ml 

n 2 



(B5) 
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FIG. 8. The left panel illustrates the temperature dependence of the potential. The continuous-red curve is tuned to give 
a vanishing vacuum energy at the minima for T = 0. Increasing the temperature shift the minima to a de Sitter vacuum 
(blue-dashed curve) and further increase of the temperature pushes the system into a symmetric high temperature phase 
(yellow-dotted curve). The right panel is illustrates the case in which the zero temperature minima is tuned to give a AdS 
vacuum (red-continuous curve), while blue-dotted curve gives a Minkowski minima and the yellow-dotted curve the de Sitter 
minima. (Left Panel) The parameters have been set to n — 100 and a = 7. The curves correspond to the following values of 
the temperature: T = 0.1 (bottom), T = 3.0 (central), T = 4.0 (top). (Right Panel) The parameters have been set to p = 120, 
a = 7. and the curves correspond to the following values of the temperature: T = 0.1 (bottom), T = 3.0 (central), T = 3.5 
(top). 



The effective potential (|B5[) has a minimum at 



Ml ~ /2 2 



1- 



1 7r 2 (8a 2 + 25) 



where 



V n 



64tt 2 



3 7r 2 (25 + 8a 2 ) 
2/I 2 a 2 



1 - 



1 7r 2 (25 + 8a 2 ) 
p? a 



2^ 1 2 



A de Sitter minimum exists in the range 



^ 2 (25 + 8a 2 ) 2 4 7r 2 (25 + 8a 2 
3 



a' d a* 

At small temperature, expanding appropriately and then summing over < I G N, we find 



eff 



25 + 8a 2 
32a 2 



Mi 
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32tt 
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2 a 
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4 ~ 




kH 2 J. 



In the high temperature limit, we find 



U e tt = VeS 



1440tt 2 
25 + 8a 2 



-Mi 



1 



+4E 



32a 2 a 32tt 2 
'x(M 2 + if 2 )' 



(=i 



7T 2 ^ 2 



| + (M 2 + i? 2 ) 2 

K-2 I 



1 fM^ + H 2 
— In 



' M 2 + H 2 



X 



(B6) 



(B7) 



(B8) 



(B9) 



(B10) 



where we have summed over I and proceeded with appropriate analytic continuation. The above expression is 
normalized by subtracting the vacuum energy contribution for M a = 0. In the adiabatic approximation adopted, the 
effect of increasing the temperature is to uplift the minimum of the potential without changing its shape around the 
minimum. 

I 
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